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(1)  to  gain  insight  into  the  fracture  behaviour  of  long  bones  under 
singular  and  combined  states  of  loading  through  a  spectrum  of  load¬ 
ing  rates,  and  we  to  characterize  the  geometric  and  mechanical 
properties  of  long  bones  by  non-destructive  means. 

Using  a  modified  Instron  Materials  Test  System  fresh 
paired  canine  long  bones  are  being  studied  under  torsional  and 
torsion-compressive  states  of  loading  through  a  spectrum  of  loading1 
rates.  Fracture  behaviour  is  characterized  by  the  torque  and 
energy  to  failure,  angle  of  twist,  fracture  type  and  location. 

Additional  studies  involve  the  development  of  techniq^A?s 
employing  computerized  axial  tomography  (CT-scan)  and  automated 
data  reduction  to  obtain  long  bone  geometric  and  mechanical  proper¬ 
ties.  Computerized  axial  tomography  allows  construction  of  sequen¬ 
tial  radiographs  in  a  plane  perpendicular  to  that  of  a  normal  flat 
plate  x-ray.  This  capability  allows  accurate  construction  of  the 
cross  sectional  dimensions  and  localized  tissue  densities  of  long 
bones  under  in  vivo  and  in  vitro  conditions.  A  correlation  has 
been  shown  to  exist  between  bone  density  and  mechanical  properties. 
To  obtain  the  geometric  properties,  a  program  called  PROPS  has  been 
developed  in  the  laboratory  and  used  in  conjunction  with  a  mini¬ 
computer  and  digitizer  to  analyze  the  scans. 

Using  these  procedures,  paired  cadaver  human  and  baboon 
femora,  tibiae,  fibulae,  humeri  and  radii  have  been  analyzed. 
Experimentally  produced  fractures  have  been  related  to  minima  of 
geometric  and  physical  properties  and  compared  with  clinically 
observed  fractures. 
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I.  INTRODUCTION 


Fractures  of  long  bones  are  a  common  clinical  finding  in 
human  trauma  situations.  These  include  injuries  which  are  sus¬ 
tained  during  automotive  and  aircraft  crashes,  parachute  land¬ 
ings  and  aircraft  ejection.  Many  of  the  latter  are  wind  flail 
induced  upper  extremity  injuries  resulting  from  hyperextension 
of  the  elbow  and  h'- oerabduction  of  the  shoulder.  Lower  extrem¬ 
ity  injuries  of  the  tibia,  fibula  and  femur  are  also  common 
findings  (Fryer,  1961;  Payne  and  Hawker,  1974;  Combs,  1977).  In 
these  situations  severe  torques,  moments  and  compressive  forces 
act  in  concert  to  produce  complex  states  of  loading.  With  the 
exception  of  work  by  Schneck  (1976,  1979),  little  information 
has  been  reported  that  attempts  to  define  the  precise  nature  of 
this  environment.  In  general,  however,  fracture  behaviour  is 
related  to  the  loading  state  and  the  rates  of  application  of  in¬ 
jury  forces,  as  well  as  local  geometrical  and  physical  proper¬ 
ties.  The  objectives  of  this  research  effort  are  twofol'd :  (1) 

to  gain  insight  into  the  fracture  behaviour  of  long  bones  under 
singular  and  combined  states  of  loading  through  a  spectrum  of 
loading  rates  in  a  controlled  laboratory  environment,  and  (2)  to 
characterize  the  geometric  and  physical  properties  of  long  bones 
by  non-destructive  means.  The  methods  developed  to  determine 
the  latter  employ  computerized  axial  tomography  and  automated 
data  reduction  and  have  resulted  in  a  significant  capability  to 
measure  the  anthropometric  properties  of  long  bone. 


II 


EFFECT  OF  LOADING  RATE  ON  FRACTURE  BEHAVIOUR 


Long  bone  fractures  usually  occur  due  to  impact  of  ex¬ 
tremities  with  stationary  objects  or  from  rapidly  applied 
loads  which  cannot  be  accommodated  due  to  soft  and  hard  tis¬ 
sue  restraints.  Since  bone  is  a  viscoelastic  material  its 
response  to  these  fracture  producing  environments  will  be  a 
function  of  loading  rate.  This  aspect  of  bone  behavior  is 
examined  in  our  study  by : 

1.  Determination  of  the  effect  of  loading  rate  upon 
bone  failure  under  torsional  loading  conditions; 

2.  Determination  of  the  effect  of  axial  compressive 
loads  on  torsional  fractures; 

3.  Analysis  of  the  effect  of  cancellous  bone  on  tor¬ 
sional  fracture  characteristics. 

2.1  Previous  Work 

Although  a  considerable  number  of  studies  have  investi¬ 
gated  the  mechanical  properties  of  standardized  specimens  and 
of  whole  bone,  most  have  been  conducted  at  low  loading  rates. 
The  following  is  a  chronological  review  of  pertinent  studies 
that  have  included  testing  at  higher  loading  rates. 

McElhaney  and  Byars  (1965),  using  standardized  speci¬ 
mens  of  bovine  and  human  femora,  conducted  constant-velocity 
compression  tests  at  strain  rates  varying  from  .001  to  1500  sec 
Plotting  energy  versus  strain  rate,  they  noted  an  increase  in 
the  energy  absorbing  capacity  of  bone,  a  peak,  then  a  sharp 
decline  as  strain  rate  increased.  They  concluded  that  the 


mechanical  properties  of  bone  are  rate  sensitive  and  that  a 
critical  velocity  for  bone  occurs  at  a  strain  rate  of  approxi¬ 
mately  1  sec”^. 

Experimental  impact  torsion  tests  of  intact  fresh  and  em¬ 
balmed  human  tibia  have  been  performed  by  Burstein  and  Frankel 
(1968).  The  tibia  was  rigidly  fastened  to  a  pendulum  and  the 
unit  was  allowed  to  gain  kinetic  energy  at  the  expense  of  its 
potential  energy.  When  maximum  velocity  was  attained,  the  grip 
on  the  distal  end  of  the  tibia  struck  a  stop  and  the  pendulum, 
with  the  proximal  end  of  the  bone,  continued  to  rotate.  The 
loading  duration  was  .1  second.  A  second  series  of  tests  were 
then  conducted  at  a  loading  duration  of  one  minute  to  determine 
if  there  were  significant  differences  in  the  energy  absorbing 
capacity.  The  tests  showed  a  considerable  variation  among  the 
individual  bones.  However,  they  were  able  to  obtain  three 
paired  specimens  thus  allowing  direct  comparisons.  A  reduction 
in  energy  capacity  of  about  55%  was  found,  compared  to  only  a 
20%  difference  noted  by  McElhaney  and  Byars  in  their  standard¬ 
ized  specimens.  They  attributed  this  difference  to  the  presence 
of  cancellous  structure  and  a  more  complex  stress  field  existing 
in  whole  bones. 

Piekarski  (1970),  using  machined  bovine  femur  specimens 
in  bending  tests,  measured  the  energy  required  for  crack 
propagation  and  observed  microscopically  the  mechanism  of 
fracture.  He  noted  that  at  low  strain  rates  cracks  propagate 
along  Haversian  lamellae  and  at  high  rates  propagate 
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indiscriminately  through  all  the  micro-constituents.  He  con¬ 
cluded  that  bone  behaves  as  a  tough  material  at  low  strain 
rates  and  fractures  as  a  brittle  material  at  high  strain 
rates . 

Sammarco ,  Burstein,  and  Frankel  (1971)  conducted  a 
study  to  determine  the  effect  of  loading  rate  on  the  ultimate 
properties  of  bone.  Their  results  were  obtained  from  tor¬ 
sional  tests  of  31  pairs  of  fresh  canine  bone  tested  at  both 
a  high  and  a  low  loading  rate.  The  right  bones  were  broken 
rapidly,  with  a  loading  time  of  40-60  milliseconds  for  frac¬ 
ture,  while  the  left  bones  were  broken  slowly,  with  loading 
times  of  between  120-180  seconds  required  for  fracture.  For 
each  pair,  and  for  each  parameter  studied  they  calculated  the 
ratio  ( F-S) / ( ( F  +  S)/2)  where  F  is  the  value  for  the  bone  loaded 
quickly,  and  S  for  that  loaded  slowly.  The  parameters  for  which 
F  and  S  could  stand  were  maximum  torque,  energy  absorption,  and 
angular  deformation.  Expressed  as  a  percentage,  angle  of  twist 
was  found  to  be  112  per  cent;  maximum  torque,  133  per  cent;  and 
energy  recorded  to  fracture,  144  per  cent.  They  concluded  that 
whole  dog  bone  exhibits  viscoelastic  properties  and  the  amount 
of  torque  and  energy  the  bone  absorbs  increased  at  the  higher 
loading  rate. 

Panjabi,  White,  and  Southwick  (1973)  conducted  a  simi¬ 
lar  torsion  experiment  on  fresh  rabbit  femora  and  ulnae.  The 
same  mechanical  parameters  were  determined  under  five  loading 
rates  ranging  from  .003  to  13.2  radians  per  second.  They  found 


that  maximum  torque,  maximum  torsional  deformation  and  energy 
absorbed  to  fracture  all  increased  with  the  rate  of  deformation, 
reached  a  maximum,  and  then  declined. 

Using  uniaxial  tensile  tests,  Wright  and  Hayes  (1976) 

studied  the  relation  between  strain  rate  and  density  on 

ultimate  strength.  They  tested  standardized  longitudinal 

specimens  of  bovine  compact  bone  at  strain  rates  ranging  from 
-4  -1 

5.3  x  10  to  237  sec  .  The  test  results  demonstrated  a 
highly  significant  correlation  between  ultimate  strength  and 
strain  rate,  and  between  ultimate  strength  and  density.  Strain 
rate  also  had  a  significant  effect  on  the  energy  absorbing 
capacity  of  bone.  The  specimens  tested  at  higher  strain  rates 
generally  absorbed  more  energy  during  failure.  However,  unlike 
some  previous  studies,  they  found  no  evidence  of  a  critical  rate 
at  which  the  energy  absorping  capacity  attains  a  maximum. 

From  a  slightly  different  viewpoint,  Ducheyne ,  Martens,  et 
al.  (1977),  investigated  the  microscopic  fracture  aspects  of  im¬ 
pact  tested  bones.  Human  tibiae  and  femora  were  torsion  tested 
to  yield  fractures  within  .05  seconds  and  30  seconds.  A  scanning 
electron  microscope  was  used  to  correlate  fracture  behavior  with 
bone  structure.  Throughout  the  impact  fracture  surfaces  both 
ductile  and  brittle  fracture  characteristics  were  observed. 
However,  the  brittle  fracture  mode  was  predominant  during  propa¬ 
gation,  indicating  a  higher  energy  absorption  in  this  fracture 
regime.  They  remarked  that,  "The  conclusion  that  a  higher 
energy  absorption  to  failure  upon  increasing  strain  rate  is 
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linked  to  an  increased  brittleness  and  a  more  marked  brittle 


fracture  surface  morphology  is  as  yet  not  established." 

A  recently  published  study  conducted  by  Martens,  et  al . 
(1980),  utilized  human  femoral  and  tib’al  whole  bone  specimens. 
These  specimens  were  subjected  to  torsional  loading  at  high 
strain  rates;  fractures  occurring  within  100  milliseconds. 

Torque  versus  angular  deflection  was  recorded.  They  calculated 
strain  energy  before  fracture,  torsional  stiffness,  maximum 
torque  and  ultimate  deflection.  A  high  variability  of  these 
mechanical  parameters  was  observed,  and  was  attributed  to 
differences  in  geometry  of  the  bones. 

One  final  relevant  study,  by  Miller  and  Piotrowski  (1974), 
was  concerned  with  the  variability  of  the  torsional  strength  of 
paired  bones.  Monkey,  dog,  and  rat  bones  were  utilized  to 
determine  if  a  right-left  bias  existed.  No  statistical  differ¬ 
ence  v/as  noted  and  it  was  concluded  that  the  torsional  strength 
of  matched  bone  pairs  can  be  expected  to  be  equal. 

2.2  Methods  / 

2.2.1  Torsional  Loading 

An  examination  of  the  effects  of  loading  rate  on  bone  un¬ 
der  a  torsional  load  has  been  conducted  in  an  attempt  to  clarify 
this  relationship.  Torsional  loading  was  chosen  for  several 
reasons:  1)  in  every  cross  section  of  a  test  body  the  torque 

remains  constant,  therefore  causing  failure  to  occur  in  the 
weakest  section;  2)  torque  and  angular  velocity  can  be  easily 
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measured  and  3)  torsion  is  a  realistic  loading  mode  producing 
failure  similar  to  clinically  observed  fractures. 

Paired  canine  radii  were  tested  utilizing  the  entire  load¬ 
ing  rate  spectrum  of  an  Instron  Testing  Mac’ ine  (Model  1115). 
This  system  is  capable  of  rates  ranging  from  .0533  to  53.3 
revolutions  per  minute  (.0558  to  5.58  radians  per  second). 

The  bones  obtained  for  this  experiment  were  from  adult 
canine  mongrels  utilized  for  various  research  projects  not 
affecting  the  skeletal  system.  As  soon  as  possible  after  ani¬ 
mal  sacrifice,  the  bones  were  removed  and  cleared  of  all  soft 
tissue.  The  bones  were  paired  and  kept  wet  until  tested.  If 
the  specimens  could  not  be  tested  the  day  they  were  obtained, 
they  were  packaged  and  frozen  at  -20°  C.  Freezing  has  been 
shown  to  have  no  significant  detrimental  effect  on  the  mechani¬ 
cal  properties  of  bone  (Sedlin,  1966),  an  observation  confirmed 
by  our  results. 

In  preparation  for  testing  the  bones  were  mounted  in  brass 
fixture  cups  with  center  pins.  Proper  position  and  alignment 
were  maintained  by  placing  the  cups  in  a  drill  press,  using  it 
as  a  vise.  Cerrobend,  an  alloy  with  a  melting  temperature  of 
158°  F,  was  then  poured  around  the  bone  ends  (Figure  1).  After 
the  Cerrobend  had  set,  the  brass  fixture  cups  were  placed  within 
a  second  set  of  concentric  cups  and  mounted  in  the  Instron. 

This  outer  pair  of  cups  was  designed  to  insure  that  the  torsion 
unit  of  the  Instron  has  accelerated  to  the  desired  rate  before 
actual  loading  of  the  bone  took  place.  This  acceleration 
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occured  as  the  outer  cup  spun  around  the  inner  cup.  The  two 
cups  engage  by  means  of  a  pin  and  the  bone  was  then  externally 


The  right  bone  of  each  pair  was  tested  at  the  slowest 
loading  rate  (.0533  rpm) ,  while  the  left  one  was  tested  at  one 
of  the  faster  rates  (.267,  .533,  2.67,  5.33,  10.33,  26.7,  53.3 

rpm).  This  pattern  was  not  altered  as  it  has  been  shown  that 
there  is  no  significant  side-to-side  variation  in  the  mechani¬ 
cal  properties  of  paired  bones  (Miller  and  Piotrowski ,  1974). 
Since  the  loading  rate  of  the  right  bone  was  kept  constant,  it 
was  used  as  a  control  to  compare  the  effect  of  increased  load¬ 
ing  rate  on  the  left  bone.  The  parameters  recorded  were:  max¬ 
imum  torque,  energy  absorbed,  angle  of  twist  at  failure,  and 
the  time  required  for  fracture.  Also  noted  was  the  type  and 
location  of  the  fracture  and,  when  possible,  the  fracture  angle 
was  measured  with  a  goniometer. 

The  Instron  load  cell  was  connected  to  a  Gould  strip 
chart  recorder  (Model  2400)  with  the  vertical  pen  displacement 
indicating  the  magnitude  of  the  torque  and  horizontal  displace¬ 
ment  indicating  the  angle  of  rotation.  Using  a  Hewlett-Packard 
9874A  digitizer  and  9845T  desktop  computer,  the  area  under  the 
torque-rotation  curve  was  calculated.  This  area  represents  the 
amount  of  energy  absorbed  by  the  bone  until  fracture.  A 
Tektronix  oscilloscope  was  also  used  to  insure  that  the  pen 
response  of  the  chart  recorder  was  sufficient  to  capture  the 
event . 
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2.2.2  Axial  and  Torsional  Loading 

A  similar  study  was  conducted  to  determine  the  effect  of 
combined  loads  at  different  loading  rates.  The  mounting  and 
testing  procedure  previously  described  in  Section  2.2.1  was 
followed  except  that  the  bone  was  subjected  to  a  compressive 
load  in  addition  to  the  applied  torque.  The  right  bone  was 
used  as  a  control  and  was  broken  at  the  lowest  loading  rate  in 
pure  torsion.  This  method  made  it  possible  to  note  the  effect 
of  the  axial  load  and  also  to  compare  these  results  with  those 
from  the  previous  study  using  torsional  loading  only.  A  com¬ 
pressive  load  of  45  kg  was  applied  axially  to  the  bones  by  a 
loading  fixture  comprised  of  three  pneumatic  cylinders  (Figure 
3).  The  bones  were  then  torsionally  loaded  utilizing  the 
entire  loading  rate  spectrum  as  before.  The  same  parameters, 
torque,  energy  and  angle  of  twist  were  recorded. 


Figure  3 

Testing  and  recording  apparatus- 
combined  torsional  and  axial  loading  experiment 
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2.2.3  Cancellous  Bone  Effects 

To  investigate  the  effect  cancellous  bone  has  on  energy 
absorption,  canine  tibiae  were  prepared  in  a  manner  similar  to 
that  previously  described  in  Section  2.2.1,  with  the  exception 
that  one  of  the  bones  from  each  pair  was  cleared  of  all  cancel¬ 
lous  bone.  This  was  accomplished  by  drilling  a  hole  in  each 
end  of  the  bone  and  reaming  the  medullary  canal.  Flexible 
reamers  capable  of  negotiating  the  curvature  of  the  shaft  were 
used  (Berry,  et  al.,  1978).  Both  the  reamed  and  unreamed  bone 
of  each  pair  was  then  torsionally  loaded  at  one  of  three  load¬ 
ing  rates  covering  the  available  spectrum  of  the  Instron  (.533, 
5.33,  and  53.3  revolutions  per  minute).  The  maximum  torque  and 
amount  of  energy  absorbed  to  failure  was  again  measured. 

Canine  tibiae  were  chosen  for  these  tests  rather  than 
radii  as  their  larger  size  and  straighter  shaft  facilitated 
the  reaming  procedure.  Additionally,  since  the  tibia  has  a 
large  amount  of  trabecular  bone,  its  absence  might  produce  a 
more  pronounced  effect. 

2.3  Results 

2.3.1  Torsional  Loading 

Y/hile  establishing  the  experimental  protocol  70  pairs  of 
specimens  were  broken.  These  preliminary  results  are  not 
reported.  Fractured  specimens  which  showed  signs  of  immaturity 
or  any  physiologic  abnormality  are  also  not  represented  in  the 
reported  data. 
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To  compare  the  effects  of  the  loading  rate  between  differ¬ 
ent  bone  pairs,  a  non-dimensional  quantity  was  calculated 
according  to  the  following  relation,  X  =  X. /XD ,  where  the 
variable  X  stands  for  one  of  the  three  parameters  recorded. 

X  indicates  the  right  bone  which  was  always  tested  at  the 
slowest  loading  rate  and  X^  indicates  the  left  bone  which  was 
broken  at  one  of  the  higher  loading  rates.  The  three  param¬ 
eters  for  which  the  variable  X  can  stand  are  the  maximum 
torque  attained  (T),  energy  absorbed  until  fracture  (E),  and 
the  angle  of  twist  at  failure  (A).  These  dimensionless  ratios 
enable  the  comparison  of  different  paired  specimens  by  minimiz¬ 
ing  the  effects  of  bone  size,  geometry,  histological  differ¬ 
ences,  or  other  factors  that  may  exist  between  different  pairs 
of  bones.  These  non-dimensional  parameters  are  an  indication 
of  the  effect  of  the  torsional  loading  rate. 

Ten  pairs  of  canine  radii  were  broken  at  each  of  the 
seven  different  loading  rates,  a  total  of  140  specimens. 

Means  and  standard  deviations  for  the  parameters  T  and  E-  were 
plotted  versus  loading  rate.  A  best-fit  curve  analysis  was 
performed  and  the  resulting  polynominal  curve  was  drawn 
(Figures  4  and  5).  Correlation  coefficients  for  T  and  E  were 
.85  and  .73,  respectively.  Means  and  standard  deviations  for 
the  angle  of  twist  data  was  plotted;  no  trend  is  evident 
(Figure  6). 


LOADING  RATE  (rtv/nin) 

Figure  4 

Torque  versus  loading  rate  (torsional  load) 


LORDING  RATE  (p»»/i»ln) 


Figure  5 

Energy  versus  loading  rate  (torsional  load) 


2.0 

1.8 

1.6 

1.4 

^  1.2 
1.0 

0.8 

0.6 

0.4 

01  0.1  1  10  100 
LORDING  RRTE  (rov/min) 

Figure  6 

Angle  of  twist  at  fracture  versus  loading  rate  (torsional  load) 


2.3.2  Axial  and  Torsional  Loading 

To  analyze  the  results  of  the  combined  loading  tests 
the  same  procedure  as  outlined  in  Section  2.3.1  for  the  tor¬ 
sional  loading  studies  was  followed.  The  parameter  indi¬ 
cates  the  left  bone  which  was  axially  loaded.  Ten  pairs  of 
canine  radii  were  broken  at  eight  different  loading  rates 
(160  specimens).  Again  a  best-fit  curve  analysis  was  per¬ 
formed  and  the  resulting  polynomial  along  with  the  means  and 
standard  deviations  for  the  parameters  T  and  E  were  plotted 
(Figures  7  and  3)  Correlation  coefficients  for  T  and  E  were 
.78  and  .86,  respectively.  As  before,  no  trend  was  realized 
in  the  angle  of  twist  data  (Figure  9). 
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Figure  9 

Angle  of  twist  at  fracture  versus  loading  rate  (torsional  axial  load 


2.3.3  Cancellous  Bone  Effects 

A  non-dimensional  ratio  was  again  used  to  compare  the 
difference  between  the  reamed  and  the  non-reamed  bones.  In 
this  case  the  ratio  used  was  X  =  where  XR  stands  for 

the  reamed  bone  and  XNR  stands  for  the  bone  that  was  not 
cleared  of  cancellous  structure.  X  represents  either  maximum 
torque  (T)  or  energy  absorbed  (E). 

The  data  was  subjected  to  a  Student's  t  test  to  deter¬ 
mine  if  cancellous  bone  plays  a  significant  role  in  energy  ab¬ 
sorption  at  the  loading  rates  used.  The  results  indicate  that 
no  significant  difference  (p<£  .25)  exists  between  the  reamed 
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and  non-reamed  bones  (Table  I).  Based  on  the  specimens  tested, 
this  finding  implies  that  long  bones  may  be  adequately  modelled 
as  hollow  tubes  rather  than  as  composite  structures. 


TABLE  I 

CANCELLOUS  BONE  EFFECTS 
CANINE  TIBIA 

Five  Pairs  at  Each  Loading  Rate 


_ TORQUE _ ENERGY 

(rev/min) 


.533 

T  = 

.916 

E 

= 

.884 

t  = 

3.89 

t 

= 

1.335 

(P  ^ 

.0025) 

(P 

.25) 

5.33 

T  = 

.960 

E 

= 

.906 

t  = 

.689 

t 

= 

.762 

(P  ^ 

.25) 

(P 

.25) 

53.3 

T  = 

.836 

E 

= 

.750 

t  = 

4.89 

t 

= 

3.786 

(P  ^ 

.001) 

(P 

.005) 

"^reamed 

T 

non-reamed 


^reamed 

^non-reamed 


17 


2.4  Discussion 


Both  pairs  of  curves,  T  and  E  versus  loading  rate  for 
torsional  loading  and  for  combined  torsional  and  axial  loading, 
are  of  similar  nature.  On  a  semi-log  plot,  all  curves  start  at 
a  minimum,  increase  slowly,  reach  a  peak  and  then  drop  off. 
However,  when  comparing  T  for  pure  torsion  to  T  for  combined 
loading,  the  combined  loading  plots  exhibit  a  slightly  higher 
magnitude,  indicating  an  effect  of  axial  load  on  the  bone. 
Maximums  in  the  curve  occured  at  15.5  rev/min  and  14.8  rev/min 
for  T  and  ~E ,  respectively,  under  pure  torsional  loading.  For 
axial  and  torsional  loading,  the  maximums  for  T  and  E  occurred 
at  17.6  rev/min  and  11.8  rev/min,  respectively. 

The  slightly  increased  torque  to  failure  which  occurs 
when  a  concurrent  compressive  axial  load  is  applied  appears 
reasonable.  Compressive  stresses  may  act  to  mitigate  torsion 
induced  tensile  stresses.  It  would  be  interesting  to  see  if 
the  effect  of  higher  axial  loads  would  be  to  exaggerate  this 
difference,  in  the  absence  of  buckling. 

Most  bones  tested  exhibited  a  spiral  type  fracture  typi¬ 
cal  of  tubular  structures  loaded  in  torsion.  Fracture  gener¬ 
ally  initiated  at  the  proximal  end  of  the  bone.  The  loading 
rate  did  not  appear  to  influence  the  position  on  the  bone  where 
fracture  occured  as  suggested  by  Sammarco ,  et  al .  (1971).  Com¬ 
minution  occured  frequently  at  the  rapid  loading  rates,  indi¬ 
cating  the  relationship  of  fracture  to  energy  dissipation, 
since  crack  formation  is  an  energy  dissipating  mechanism 
( Figure  10 )  . 
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Figure  10 

Spiral  Fracture  (top)  and  Comminution  (bottom) 


Torsional  testing  conducted  on  a  large  population  over  a 
wide  range  of  loading  rates  is  not  found  in  the  literature. 
Several  investigators  have  noted  an  increase  in  the  energy 
absorbing  capacity  of  bone  at  higher  loading  or  strain  rates. 
Sammarco,  et  al.,  found  this  increase,  testing  specimens  at 
only  two  rates,  thereby  implying  a  continuously  increasing 
energy  absorption.  Panjabi,  et  al .  (1973),  using  five  differ¬ 
ent  loading  rates  obtained  curves  for  torque  or  energy  versus 
loading  rate  that  increased  to  a  maximum  and  then  declined. 
This  trend  is  in  general  agreement  with  the  results  obtained 
for  the  pure  torsion  and  the  combined  loading  tests  conducted 
in  our  study.  Panjabi  averaged  the  results  from  both  rabbit 
femora  and  ulnae  but  did  not  indicate  the  statistical  signifi¬ 
cance  of  the  data.  In  our  study,  1G0  specimens,  ten  pairs  at 
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each  of  eight  selected  loading  rates,  were  tested  producing 
curves  with  relatively  high  correlation  coefficients.  By  using 
a  large  sample  population  and  non-dimensional  parameters,  ex¬ 
perimental  artifacts  and  biological  variations  are  minimized. 

Studies  relevant  to  the  fracture  mechanisms  of  whole 
bones  under  combined  torsional  and  axial  loading  and  at  differ¬ 
ent  loading  rates  have  not  been  previously  reported.  Fractures 
produced  under  this  combined  state  of  loading  are  representa¬ 
tive  of  those  clinically  found  in  human  trauma  situations. 


III.  TOPOGRAPHIC  PHYSICAL  AND  CROSS  SECTIONAL  PROPERTIES 


The  fracture  behavior  of  a  long  bone  cannot  be  completely 
characterized  unless  the  geometry  of  the  bone  and  any  regional 
variations  in  strength  within  the  bone  are  determined.  An  ideal 
method  to  accomplish  this  would  be  non-destructive  so  that  the 
specimens  would  also  be  available  for  histological  study  or  me¬ 
chanical  testing.  The  data  derived  from  such  experiments  could 
be  put  to  other  uses,  such  as  modelling  of  bone  fractures. 

In  this  section  the  geometric  and  physical  properties  of 
long  bones  have  been  determined.  One  of  the  exciting  fallouts 
of  this  work  has  been  the  development  of  non-destructive  tech¬ 
niques  employing  computerized  axial  tomography  (CT-scan)  and 
automated  data  reduction  to  obtain  this  information.  Originally 
it  had  been  proposed  to  utilize  hardness  plots  to  determine 
local  strength  variations.  This  method  can  only  be  employed  on 
the  bone  surface  and  due  to  the  curvature  of  the  bone  repeata¬ 
bility  is  poor.  Additionally,  the  Structural  Dynamics  Research 
Corporation  SPOTS  program  was  originally  used  to  determine  the 
cross  sectional  properties.  It  was  found  to  be  time  consuming 
and  tedious  in  application,  and  to  have  difficulty  representing 
certain  bone  cross  sections.  Computerized  axial  tomography 
allows  construction  of  sequential  radiographs  in  a  plane  perpen¬ 
dicular  to  that  of  a  normal  flat  plate  x-ray.  This  capability 
permits  accurate  simultaneous  construction  of  the  cross  sec¬ 
tional  dimensions  and  delineation  of  localized  tissue  densities 
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of  long  bones  under  in  vivo  and  in  vitro  conditions.  In  addi¬ 
tion,  the  data  is  stored  on  a  magnetic  medium  and  may  be  readily 
recalled  and  manipulated.  A  method  has  been  developed  to 
rapidly  and  accurately  process  the  vast  amount  of  data  which  is 
generated . 

3.1  Previous  Work 

Of  the  seven  articles  found  in  the  literature  on  prop¬ 
erty  determination  of  bone,  three  are  concerned  with  the  geo¬ 
metric  properties  of  the  human  tibia,  one  deals  with  physical 
properties  of  the  human  femur,  while  the  remaining  three  fur¬ 
ther  develop  algorithms  and  procedures  to  compute  cross  sec¬ 
tional  geometric  properties. 

The  first  of  the  papers  (Minns,  Bremble  and  Campbell, 

1975)  described  a  method  in  which  seven  serial  sections  are 
taken  at  intervals  ranging  from  20  to  80%  of  the  total  bone 
length.  A  cartesian  grid  is  superimposed  upon  a  photograph 
of  the  section  and  the  coordinates  of  the  grid  ai'e  used  as  in¬ 
put  for  a  computer  program  to  calculate  the  geometric  proper¬ 
ties.  The  results  from  four  left  human  tibiae  showed  that  the 
minimum  cross  sectional  area  occurred  at  80%  of  the  total 
length.  The  minimum  moment  of  inertia  fell  between  70  and  80% 
of  the  length  and  the  minimum  polar  moment  of  inertia  occurred 
at  70%  where  it  would  be  expected  that  mechanical  failure  due 
to  torsional  loading  only  would  occur. 
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The  second  of  the  papers  on  the  geometric  properties  of 
the  human  tibia  (Miller  and  Purkey,  1980)  analyzed  right-left 
variability,  sex  related  differences,  and  expected  difference 
between  individuals.  Nine  pairs  of  embalmed  human  tibiae  were 
obtained:  three  male,  two  female,  and  four  of  unknown  gender. 

Fifteen  serial  sections  were  cut  from  each  tibia  representing 
cross  sections  from  20  to  80%  of  the  total  length.  A  dimension 
reference  was  marked  on  each  section  and  an  enlarged  tracing  of 
the  cortical  bone  and  reference  mark  was  produced  with  the  aid 
of  an  opaque  projector.  The  coordinates  of  the  inside  and  out¬ 
side  boundary  of  each  cross  section  were  used  as  input  for  a 
computer  program.  The  results  indicated  a  significant  differ¬ 
ence  in  properties  between  pairs  of  bones  from  different  indi¬ 
viduals  yet  the  graphs  had  similar  shapes.  No  statistically 
significant  right-left  bias  was  observed.  Also,  the  minimum 
polar  moment  of  inertia  occurred  at  70%  of  the  length.  This 
corresponds  to  the  primary  location  of  clinically  observed 
torsional  fractures. 

Incidental  to  their  finite  element  analysis  of  the  human 
tibia,  Piziali  and  Hight  (1976)  determined  the  cross  sectional 
properties  of  four  cadaver  specimens.  Inertial  properties  were 
largest  at  the  metaphyses  and  diminished  in  the  shaft.  Cross 
sectional  area  followed  a  similar  pattern  but  the  overall  vari¬ 
ation  was  not  as  pronounced.  Minima  of  inertia  and  area 
occurred  at  approximately  20%  superior  to  the  distal  end. 
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In  a  study  of  the  physical  properties  of  the  human  femur, 
Evans  and  Lebow  (1953)  used  samples  of  compact  bone  which  were 
cut  longitudinally  from  anterior,  posterior,  medial,  and  later¬ 
al  quadrants  of  the  proximal,  middle,  and  distal  thirds  of  the 
femoral  shaft.  Both  wet  and  dry  samples  were  used,  each  sample 
being  machined  to  a  standard  size.  These  samples  were  then 
tested  to  determine  the  modulus  of  elasticity,  ultimate  tensile 
strength  and  hardness.  Femurs  from  five  male  Caucasian  cadav¬ 
ers  ranging  in  age  from  58  to  81  were  studied.  The  results  in¬ 
dicated  no  significant  difference  between  modulus  of  elasticity 
and  hardness  for  the  wet  and  dry  samples  from  various  quad¬ 
rants.  Further,  the  middle  third  exhibited  the  highest  average 
modulus  of  elasticity  (14.6  GPa),  maximum  hardness,  and  highest 
average  ultimate  tensile  strength  (83.2  MPa).  The  proximal 
third  had  the  lowest  average  modulus  of  elasticity  (13.5  GPa), 
minimum  hardness,  and  the  lowest  average  ultimate  tensile 
strength  (77.6  MPa).  Air  drying  the  samples  increased  the 
modulus  of  elasticity  17.6%  and  the  hardness  54.3%.  The 
average  ultimate  tensile  strength  also  increased  by  31.3%. 

The  final  three  papers  deal  with  development  of  geometric 
property  determination  procedures.  In  the  first  of  these 
articles,  Lovejoy  and  Burstein  (1977)  obtained  laminograins  of 
five  left  human  tibiae,  taken  at  half  their  length.  These  sec¬ 
tions  were  removed  from  the  bone  and  photographs  were  projected 
with  laminograms  onto  a  cartesian  grid.  The  coordinates  from 
the  grid  were  used  as  input  for  a  computer  program  to  calculate 


I 


24 


the  various  geometric  properties.  Laminogram  and  bone  section 
results  were  then  compared  and  the  percent  error  determined. 

The  average  error  for  all  geometric  properties  (i.e.  area, 
moments  of  inertia,  and  polar  moment)  were  less  than  6%.  This 
paper  confirms  the  accuracy  of  using  x-ray  cross  sections  of 
bone  for  geometric  property  computations.  Also,  the  procedure 
using  laminograms  allows  the  bone  specimen  to  be  left  intact  for 
further  testing. 

A  similar  study  was  done  by  the  same  author  (Lovejoy  and 
Barton,  1980)  using  a  radial  grid  instead  of  the  cartesian  grid 
for  data  reduction.  The  average  error  was  found  to  be  less  than 
8%  for  the  various  properties.  The  benefit  of  this  method  was 
that  individual  sectional  properties  could  be  determined  more 
rapidly  than  with  the  previous  cartesian  grid  method.  The 
average  time  for  data  reduction  was  reduced  from  more  than  45 
minutes  to  about  two  minutes.  With  the  advent  of  automated 
techniques,  however,  this  method  must  still  be  considered  slow. 

Nagurka  and  Hayes  (1980)  developed  an  algorithm  to 
determine  geometric  properties  of  irregular  shapes  and  suggested 
its  application  to  long  bone  studies. 

It  is  anticipated  that  knowledge  of  the  locations  of  mate¬ 
rial  and  cross  sectional  weakness,  when  correlated  to  the  loca¬ 
tion  of  fracture  sites,  will  increase  our  understanding  of  long 
bone  fracture  mechanics.  In  turn,  this  may  lead  to  better  meth¬ 
ods  of  fracture  prevention. 
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3.2  Methods 


3.2.1  Cross  Sectional  Property  Software 

An  interactive  program  has  been  written  to  determine  the 
bone  cross  sectional  properties.  The  complexity  of  the  shapes 
involved  precludes  manual  computation  of  these  properties. 
Originally,  as  has  been  mentioned,  it  was  intended  to  use  SPOTS 
(Section  Properties  of  Thin-walled  Sections)  to  facilitate  sec¬ 
tional  property  data  collection.  However,  the  amount  of  manual 
pre-processing  necessary  to  faithfully  represent  the  sections 
was  found  to  be  considerable.  The  SPOTS  program  is  finite  ele¬ 
ment  based  and  requires  the  bone  to  be  broken  up  graphically  in¬ 
to  rectangles  and  arc  segments.  Certain  coordinates  of  these 
base  figures  are  then  used  to  define  the  bone  shape.  The  con¬ 
struction  of  many  small  elements  is  necessary  for  the  complex 
shapes  involved. 

Because  of  the  tedious  nature  of  the  above  procedure,  it 
was  decided  to  approach  the  problem  from  a  different  viewpoint. 

A  Hewlett-Packard  desktop  computer  (Model  9845T)  and  digitizer 
(Model  9874A)  were  used  to  implement  section  property  software 
based  partially  on  an  algorithm  presented  by  Wojciechowski 
(1976).  The  fundamental  equations  and  a  basic  flow  chart  of 
this  program  called  PROPS  are  contained  in  Appendix  I .  The  sec¬ 
tion  under  consideration  must  be  entirely  contained  in  the  first 
quadrant  of  the  digitizer  coordinate  system.  Coordinates  may  be 
digitized  at  intervals  as  small  as  .1  mm,  but  a  1  mm  interval 
has  been  found  to  give  reliable  results. 


The  procedure  is  fast 


and  efficient  since  the  operator  merely  traces  the  exterior  bound 
ary  of  the  bone  in  a  clockwise  direction  and  the  interior  medul¬ 
lary  canal  in  a  counterclockwise  direction  with  the  digitizer 
cursor.  Pertinent  cross  sectional  property  output  includes  area, 
maximum  and  minimum  principal  moments  of  inertia  and  polar  moment 
of  inertia. 

The  analytical  solutions  for  properties  of  selected  shapes 
have  been  compared  with  results  numerically  calculated  by  PROPS. 
It  can  be  seen  by  the  comparison  in  Table  II  that  agreement  is 
excellent . 


TABLE  II 


Percent  error  in  selected  properties. 


The  only  anomaly  that  has  been  found  is  that  when  two  successive 
points  have  the  same  x  coordinate,  a  division  by  zero  occurs  in 
the  product  of  inertia  calculation.  The  program  notes  this  occur 
rence ,  ignores  the  term  in  the  summation  and  continues  execution. 

The  PROPS  program  as  we  have  implemented  it  has  several  ad¬ 
vantages  over  SPOTS.  It  does  not  require  any  manual  measurements 
therefore  reducing  the  occurrence  of  operator  error  and  overall 
computation  time.  Much  more  data  is  acquired  and  the  nature  of 
the  algorithm  results  in  a  more  realistic  model  of  the  cross  sec- 


3.2.2  CT  Scanrun;’  Procoduro 

Spec  incus  are  kept  in  a  freezer  at  -20  °C  after  collection. 
The  desired  specimen  is  removed  from  the  freezer  to  thaw  before 
testing*  It  is  then  stripped  of  flesh  and  bony  landmarks  which 
denote  the*  area  to  be  scanned  are  identified.  The  bone  is 
sealed  in  a  circular  tank  completely  filled  with  water.  The 
water,  in  addition  to  preventing  dehydration,  acts  as  one  of  the 
references  needed  to  relate  CT  number  to  density.  The  bone  is 
then  scanned  in  10  m  increments  utilizing  an  Ohio  Nuclear  Model 
2020  scanner. 

Figure  11a  shows  a  scan  of  paired  baboon  humeri  (top)  and 
femora  (bottom).  The  scans  are  also  displayed  simultaneously  in 
color  (Figure  lib).  The  circular  objects  in  the  center  of  the 
picture  are  phantoms  used  to  check  machine  calibration  and  for 
reference  dimensions.  The  display  has  been  processed  so  that 


Color  and  black  and  white  scans  of  paired  baboon  humeri  (top)  and 
femora  (bottom) .  The  black  and  white  scan  is  5  cm  more  proximal. 
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CT  numbers  above  450  appear  white  and  those  below  450  uppeu  r 
black.  This  sharpens  resolution  of  the  boundaries  and  aids  dig¬ 
itization.  The  boundaries  are  non-distinet  due  to  partial  vol¬ 
ume  artifact,  that  is,  the  picture  elements  (pixels)  on  the 
boundary  not  only  encompass  bone  but  surrounding  water  as  well. 


The  result  is  an  artificially  low  CT  number  somewhere  between 
that  of  bone  and  water.  This  appears  as  a  light  blue  halo  in 
Figures  ll-14b.  Since  the  minimum  CT  number  of  cortical  bone  is 
approximately  900  and  that  of  water  is  0 ,  a  cutoff  value  of  450 
rejects  pixels  which  are  less  than  50)j  bone.  It  is  obvious  that 
the  boundaries  in  Figures  ll-14b  are  much  less  distinct  than 
those  of  Figures  ll-14a. 

A  bone  is  selected  by  positioning  a  frame  over  it 
(Figure  lib)  and  is  then  enlarged  (Figure  12).  The  coarse 
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Erroi’  can  be  introduced  if  the  bone  is  not  scanned  perpen¬ 
dicular  to  its  long  axis.  Every  effort  is  made  to  insure  that 
this  error  does  not  occur  but  a  certain  amount  is  inescapable 
since  most  long  bones  have  a  slight  natural  curvature.  An  attempt 
was  made  to  quantify  this  error.  If  a  long  bone  is  assumed  to  be 
a  thick  walled  circular  cylinder  and  if  it  is  cut  or  scanned  at 
some  angle  Q  to  its  long  axis,  a  thick  walled  elliptical  section 
will  be  revealed.  The  properties  of  the  circular  and  elliptical 
sections  have  been  compared  as  a  function  of  the  angle Q  .  The 
difference  in  area  was  found  to  be  a  function  of  sec Q  and  the 
difference  in  moments  of  inertia  was  found  to  be  a  function  of 

O 

sec  £)  .  This  is  plotted  in  Figure  15.  It  can  be  seen  that  for 
angles  as  large  as  10°  the  error  is  small. 
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Figure  15 

Estimate  of  error  introduced  by  inclination  of  the  scanning  plane 

to  the  bone  axis. 
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3.2.3  Linear  Dependency  of  Density  and  Mechanical  Properties 
on  Measured  CT  Number 

Computed  tomography  (CT)  is  a  well  developed  diagnostic 
imaging  method  for  determining  the  linear  attenuation  coeffi¬ 
cient  of  x-rays  passing  through  an  object.  An  image  of  numbers 
which  is  proportional  to  the  linear  attenuation  coefficients  is 
computed  from  precise  measurements  at  many  positions  and  many 
angles.  The  linear  relationship  between  these  numbers  and  the 
linear  attenuation  coefficient  is  the  defining  equation  for  the 
so-called  Hounsfield  units.  This  is  given  by 


u  .  -  u  , 

1  water 

u 

water 


where  CI\  is  the  calculated  CT  number  in  the  it!l  pixel,  is  the 
effective  linear  attenuation  coefficient  for  that  pixel,  uwater  is 
the  effective  linear  attenuation  of  water,  and  k  is  the  propor¬ 
tionality  constant  which  defines  Hounsfield  units  (and  thus  cali¬ 
brates  the  CT  scanning  apparatus)  and  should  have  the  value  1000. 
On  the  Hounsfield  scale,  air  has  the  value  of  -1000,  water  is  0, 


and  bone  is  about  1000  and  above.  However,  the  constant  of 


proportionality  k  is  experimentally  established  by  calibration  of 
the  CT  scanning  apparatus  and  may  not  be  treated  as  a  universal 
constant.  It  is  effectively  constant  but  not  necessarily  equal  to 


1000  for  a  given  instrument. 

The  linear  attenuation  coefficient,  u,  strongly  depends  on 
the  effective  energy  of  the  x-ray  beam  and  the  elemental  composi¬ 
tion  of  the  material  in  the  pixel.  For  a  material  of  uniform 


composition,  tin:  linear  attenuation  coefficient  is  given  by 


where 


ui  "PN!  <(?R+  0'e+0'pe) 

is  the  density,  N  is  Avagadro’s  number,  Z  the  atomic  num¬ 
ber,  A  the  atomic  weight,  and  the  (J's  are  the  scattering  cross 
sections  per  electron  for  Rayleigh  scattering,  Compton  scattering, 
and  the  photoelectric  effect.  The  Compton  cross  section  is  mildly 
dependent  upon  the  effective  energy  of  the  x-ray  beam  and  indepen¬ 
dent  of  Z.  The  cross  sections  for  Rayleigh  and  photoelectric 
scattering  are  strongly  dependent  on  both  energy  and  atomic 
number. 

Determining  calibration  standards  from  materials  with  atomic 
numbers  and  densities  similar  to  bone  allows  removal  of  the  depen¬ 
dence  on  the  Hounsfield  scaling  constant  and  the  effective  linear 
attenuation  of  water  to  give  the  relationship 

P  i  =  C^  A  +  B 

The  constants  A  and  B  depend  on  the  measured  CT  numbers  of  the 
standards,  the  effective  energy  of  the  x-rays  and  on  the  ratios  of 
the  total  electronic  scattering  cross  sections  for  the  measured 
and  standard  materials.  This  equation  demonstrates  a  linear  rela¬ 
tionship  between  density  and  CT  number  in  our  region  of  interest. 

The  corresponding  relationship  between  density  and  mechani¬ 
cal  properties  is  well  documented.  Abendshein  and  Hyatt  (1970) 
found  a  linear  correlation  (r  =  .94)  between  wet  density  and  modu¬ 
lus  of  elasticity  for  52  specimens  of  human  ti’oial  and  femoral 
cortical  bone.  Martin  and  Atkinson  (1977)  also  determined  a 


33 


linear  relationship  (r  =  .99)  between  wet  density  and  ultimate 
bending  strength  for  37  specimens  of  human  femoral  bone. 

Evans  (1973)  studied  embalmed  human  tibial ,  femoral  and 
fibular  specimens.  For  70  tibial  specimens,  positive  correla¬ 
tions  at  the  .001  level  of  significance  were  found  between  dry 
density  and  ultimate  tensile  strength,  modulus  of  elasticity 
and  Rockwell  hardness.  At  the  .05  level  of  significance  a  cor¬ 
relation  was  found  between  dry  density,  single  shear  strength 
and  Rockwell  hardness.  For  17  femoral  specimens  a  correlation 
was  only  found  between  dry  density  and  Rockwell  hardness.  This 
was  attributed  to  the  relatively  small  population  studied. 

In  summary,  since  linear  relationships  exist  between  CT 
number  and  density  and  between  density  and  mechanical  proper¬ 
ties  a  linear  relationship  between  CT  number  and  mechanical 
properties  can  be  assumed.  Thus,  a  local  CT  number  minimum 
should  indicate  a  local  mechanical  property  minimum. 

3.2.4  Physical  and  Cross  Sectional  Data  Collection  and 
Reduction 

/ 

Slides  of  the  scans  are  back  projected  on  the  digitizer 
platen  to  determine  bone  cross-sectional  properties.  The  bone 
perimeter  is  digitized  and  PROPS  is  utilized  to  perform  the 
property  calculations.  The  graphics  capabilities  of  the 
Hewlett-Packard  desktop  computer  are  employed  to  plot  these 
properties  as  a  function  of  longitudinal  bone  position. 

CT  numbers  are  obtained  using  a  PDP-15  image  processing 
system  and  a  line  printer.  A  bone  is  displayed  on  the  CRT  and 
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The  numbers 


location  whose  CT  numbers  are  desired  is  selected, 
are  then  output  to  a  line  printer.  Maxinun  values  on  the  lateral, 
medial,  posterior  and  anterior  aspects  of  the  bone  are  recorded. 
These  values  are  also  plotted  as  a  function  of  longitudinal  bone 
position . 

3.3  Results 

The  long  bones  of  a  human  and  a  baboon  cadaver  have  been 
studied,  542  diaphyseal  sections  being  evaluated.  Complete  physi¬ 
cal  and  cross  sectional  properties  are  graphically  presented  in 
Appendix  II.  Unfortunately,  complete  CT  number  data  on  the  human 
bones  is  not  available  since  a  scanner  capable  of  differentiating 
between  high  bone  densities  was  obtained  only  recently.  For  ease 
of  comparison  some  selected  properties  of  human  and  baboon  humeri 
and  femora  are  shown  in  Figures  16-19.  Note  that  the  data  for  the 
baboon  is  less  irregular  than  that  of  the  human  since  it  was  ac¬ 
quired  at  a  later  date  and  employed  a  new  scanner  with  increased 
resolution  and  the  smoothing  and  filtering  techniques  described  in 
Section  3.2.2. 


3.4  Discussion 

Based  on  the  number  of  specimens  available,  the  following 
general  observations  are  made.  On  the  average,  minimum  cross  sec¬ 
tional  areas  of  human  bones  are  2.5  times  greater  than  that  of 
baboons,  the  difference  being  more  pronounced  in  the  lower  extrem¬ 
ities.  Minimum  inertial  properties  are  approximately  5  tines  as 
great  in  man,  again  the  difference  being  larger  in  the  lower 
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Figure  17 

Comparison  of  Human  and  Baboon 
Humeral  Polar  Moments  of  Inertia 
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Figure  19 

Comparison  of  Human  and  Baboon 
Femoral  Polar  Moments  of  Inertia 
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extremities.  This  is  probably  due  to  the  necessity  of  having  in¬ 
creased  bony  development  in  the  ape  to  complement  its  upper  body 
strength.  No  left-right  differences  are  apparent  for  man  or 
baboon.  In  most  instances  cross  sectional  properties  are  at  a 
minimum  in  the  mid-shaft  region  and  increase  at  the  metaphyses. 

An  attempt  was  made  to  compare  our  cross  sectional  property 
results  to  previous  work.  Although  many  researchers  have  studied 
techniques  to  determine  these  properties,  few  have  actually 
applied  their  methods  and  none  have  studied  all  of  the  skeletal 
long  bones  in  such  detail. 

Piziali,  Flight  and  Nagel  (1976)  presented  data  on  four 
human  tibiae  which  closely  matches  our  results  both  in  trend  and 
magnitude.  Miller  and  Purkey  (1980)  in  a  study  of  nine  pairs  of 
human  tibiae,  and  Minns,  Bremble  and  Campbell's  (1975)  study  of 
four  human  tibiae  presented  similar  results  which  are  in  general 
slightly  lower  in  magnitude  than  the  current  study  and  that  of 
Piziali  et  al .  This  is  probably  due  to  a  slightly  different  in¬ 
terpretation  of  the  location  of  the  cortical-cancellous  boundary. 
Since  visual  determination  of  this  boundary  is  quite  subjective, 
it  is  believed  that  our  technique  as  illustrated  in  Figures 
ll-14a  results  in  more  consistent  delineation. 

CT  numbers  are  generally  greatest  in  the  mid-shaft  region 
where  the  anterior  and  medial  portions  are  higher  than  those  pos¬ 
terior  and  lateral.  The  femur  and  tibia  have  the  highest  CT  num¬ 
bers  while  the  radius  and  fibula  have  the  lowest.  The  low  CT 
numbers  for  the  fibula  may  be  due  to  its  small  size  and  the 
resulting  partial  volume  artifact.  Little  right-left  difference 


is  apparent 
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The  non-dimensional  quantities  T  and  E  were  calculated  to  be 
1.36  and  1.96,  respectively.  T  agreed  quite  well  (3%  differ¬ 
ence)  with  the  average  value  of  1.32  for  the  canine  radii.  E 
was  higher  than  the  canine  average  of  1.40  but  much  more  scat¬ 
ter  is  experienced  in  the  energy  data.  The  difference  in  E  is 
mainly  accounted  for  by  a  large  angle  of  twist.  Like  the 
canine  radii  fracture  occured  at  the  same  location  on  the  prox¬ 
imal  end  of  the  bone  and  was  not  shifted  by  changing  the  load¬ 
ing  rate.  It  is  difficult  to  draw  a  conclusion  based  on  this 
experiment  alone  but  it  appears  that  it  may  be  possible  to  draw 
a  direct  or  scaled  correlation  between  species. 

Figure  21  shows  the  location  of  the  fracture  site  in  re¬ 
lation  to  CT  numbers  and  polar  moment  of  inertia.  It  can  be 
seen  that  the  fracture  zone  occurs  at  a  point  of  minimum  polar 
moment  of  inertia. 

Coincidentally,  the  left  baboon  fibula  showed  evidence  of 

a  fracture  in  the  process  of  healing.  This  is  demonstrated  by 

the  callus  formation  shown  on  the  baboon  fibula  cross  sectional 

/ 

property  graphs  in  Appendix  II.  Figure  22  shows  a  comparison 
of  polar  moment  of  inertia  and  the  CT  numbers  of  the  normal 
(right)  fibula.  Polar  moment  of  inertia  is  fairly  constant 
along  the  length  but  at  the  fracture  location  a  CT  number  mini¬ 
mum  exists  on  the  posterior  surface  of  the  bone.  It  appears 
that  polar  moment  of  inertia  controlled  the  fracture  site  loca¬ 
tion  on  the  radius  but  CT  number  controlled  the  location  on  the 
fibula,  assuming  that  the  primary  fibular  load  was  torsion. 
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Figure  22 

Callus  Formation  Indicating  the  Location  of  a  Fracture 
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This  implies  that  for  a  given  bone  and  loading  state  mechanical 
or  cross  sectional  properties  or  both  may  be  determinants  of 
failure . 

The  combined  effect  might  be  determined  by  means  of  a 

"Fracture  Score"  which  quantifies  the  relative  magnitude  of 

cross  sectional  and  material  strength.  For  example,  in  the 

most  basic  form,  by  assigning  a  value  of  one  to  the  maximum 

polar  moment  of  inertia  (Jraax)  and  average  maximum  CT  number 

(CT  )  and  a  value  of  zero  to  the  minima  of  these  parameters 

(J  .  ,  CT  .  ),  the  Fracture  Score  (FS)  at  any  point  x  along  a 
v  mm  ’  min '  ’  v  J  & 

bone  in  torsion  might  be  computed  by: 


FS 


J  -  J  .  + 

x  mm 


CT  -  CT  . 
x  mm 


J  -  J  . 
max  mm 


CT  -  CT  . 
max  min 


FS  would  have  a  value  of  1  if  x  corresponded  to  a  point  of  max¬ 
imum  CT  number  and  polar  moment  and  0  if  it  corresponded  to  a 
point  of  minimum  CT  number  and  polar  moment.  Using  this  scheme 
FS  was  calculated  for  the  left  baboon  radius  in  Figure  21  and 
plotted  (Figure  23).  This  indicates  the  probability  of  frac- 
ture  in  torsion  at  the  proximal  end  which  did,  in  fact,  occur. 

The  proceeding  example  was  done  for  illustration  and 
could  only  be  adopted  for  use  if  proven  by  comparison  with  much 
more  experimental  data.  Fracture  Score  may  have  to  accommodate 
different  combinations  of  variables  for  various  loading  states 
and  may  have  to  be  weighted  in  some  manner.  However,  it  is  be¬ 
lieved  the  utility  of  the  concept  has  been  demonstrated. 
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Figure  23 


Tabulation  of  a  Candidate  Fracture  Score 
for  the  Baboon  Radius  of  Figure  21 


4.2  Modelling 

4.2.1  The  Canine  Fracture  Model 

The  canine  radius  was  chosen  for  the  loading  rate  studies 
since  it  is  easily  mounted  and  a  large  number  of  fresh  paired 
bones  were  available.  The  primate  e::trernities  that  were  ob¬ 
tained  sometimes  were  not  paired  or  had  been  damaged  during 
dissection.  It  was  felt  that  since  the  tests  result  in  non-di¬ 
mensional  data  they  determine  the  behavior  of  bone  tissue  and 
the  data  might  be  correlated  to  other  species.  This  was  demon¬ 
strated  by  fracture  of  the  paired  baboon  radii  as  described  in 
Section  4.1.  The  baboon  bones  fractured  in  the  same  location 
as  the  canine  bones  and  T  was  in  close  agreement.  E  was 


considerably  larger  for  the  baboon,  but  the  energy  data  experi¬ 
ences  much  more  scatter  (Figures  4  arid  5). 


In  an  interesting  related  study,  Strongwater,  Goel , 
Panjabi,  and  Drinker  (1D80)  investigated  the  anatomical  simi¬ 
larity  of  human,  baboon,  dog,  and  sheep  femora.  They  noted  a 
greater  external  similarity  between  dogs  and  humans  and  greater 
internal  similarity  between  baboons  and  humans.  They  concluded 
that  either  dogs  or  baboons  are  suitable  implant  models,  but 
due  to  cost,  availability,  and  ease  of  handling,  the  dog  may  be 
more  attractive. 

By  their  nature,  biological  structures  are  considerably 
more  variable  than  engineering  structures,  both  in  geometry  and 
composition.  This  necessitates  the  evaluation  of  large  sample 
populations  to  draw  statistically  significant  conclusions. 
Specimen  availability  is  of  primary  importance. 

4.2.2  l.' t  i  1  ir.a  t  ion  of  Computed  Tomography  in  Modelling 

One  of  the  more,*  obvious  uses  of  the  computed  tomography 
derived  data  is  in  the  development  of  long  bone  models.  These 
models  may  be  used  for  long  bone  fracture  studies  or  in  the  de¬ 
sign  of  orthopaedic  appliances.  More  advanced  models  would  un¬ 
doubtedly  employ  the  finite  clement  method .  The  techniques 
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might  be  the  determination  of  body  segment  distributed  mass 
properties.  If  the  volume  (slice  thickness  x  pixel  area)  and 
density  of  a  pixel  are  known,  then  the  mass  of  the  pixel  can  be 
computed.  The  sum  of  all  these  masses  would  equal  the  mass  of 
the  whole  object.  This  would  not  require  a  major  programming 
effort  since  the  scan  information  may  be  thought  of  as  stored 
in  a  3  dimensional  array  where  each  element  is  the  CT  number  of 
a  known  pixel  volume  within  the  object.  When  a  direct  conver¬ 
sion  from  CT  number  to  density  has  been  derived  it  is  intended 
to  test  this  method  by  predicting  the  weight  of  a  bone  which 
has  been  scanned. 

4.3  Relation  to  Windflail  Injury 

Current  knowledge  of  the  loads  to  which  the  human  musculo¬ 
skeletal  system  is  exposed  in  normal  activity  is  limited. 
State-of-the-art  techniques  to  obtain  this  information  generally 
combine  a  controlled  laboratory  simulation  f  the  activity  and 
theoretical  analysis.  Results  vary  among  investigators  even  un¬ 
der  these  conditions.  The  difficulties  involved  in  determining 
the  loading  state  acting  on  a  human  body  during  ejection  from  an 
aircraft  car:  therefore  be  appreciated. 

The  first  ejection  injury  occurred  over  36  years  ago 
(Fleming,  1979).  However,  due  to  its  complex  nature,  the  load¬ 
ing  environment  remains  largely  undefined.  Schneck  (1979)  cal¬ 
culated  dislodging  force  for  a  forearm  on  a  restraining  surface 
at  an  altitude  of  10,000  feet  and  various  angles  of  attack.  He 
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determined  that  dislodgment  would  occur  at  approximately  450 


knots  since  the  average  grip  retention  strength  of  270  kg  is 
exceeded.  This  agrees  well  with  Combs'  study  (1977)  of  severe 
(windflail)  and  minor  ejection  injuries.  Minor  injuries 
occurred  at  an  average  of  310  knots  while  the  windflail  inju¬ 
ries  occurred  at  an  average  of  403  knots.  Two  mid-shaft  com¬ 
minuted  fractures  of  the  femur  were  reported.  The  cross 
sectional  property  data  (Appendix  II)  for  the  human  femur  is  at 
a  minimum  between  50-60°)  of  the  scanned  length.  Additionally, 
five  proximal  humerus  fractures  were  reported.  The  cross  sec¬ 
tional  properties  determined  in  our  study  remained  constant 
along  the  length,  but  density  data,  when  available,  might 
reveal  the  reason  for  the  proximal  fractures. 

Comparison  of  fractures  produced  under  known  loading  en¬ 
vironments  to  those  produced  during  ejection  may  be  the  most 
practical  method  to  determine  a  lower  bound  on  ejection  forces 
and  moments. 
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V 


2) 


3) 


4) 


5) 


6) 


7) 


V.  SUMMARY 

The  following  outlines  the  primary  results  of  studies  performed 
in  this  contract. 

1)  For  canine  radii  in  pure  torsion,  increased  loading 
rate  produced  increased  torque  and  energy  to 
fracture.  However,  at  the  high  end  of  our  available 
loading  regime  a  drop  in  these  parameters  was  noted. 

For  canine  radii  loaded  in  torsion  and  compression 
similar  behavior  was  experienced  with  the  exception 
that  torque  and  energy  to  fracture  had  slightly 
larger  magnitudes. 

The  removal  of  cancellous  bone  was  found  not  to  have 
as  statistically  significant  effect  on  torque  and 
energy  to  fractue  or  an  effect  on  fracture  site 
location . 

A  technique  to  non-destructively  evaluate 
cross-sectional  and  physical  properties  of  long  bones 
has  been  developed. 

This  technique  has  been  used  to  analyze  the  long 
bones  of  a  human  and  a  baboon  cadaver. 

Physical  and  cross  sectional  properties  for  a  pair  of 
baboon  radii  have  been  related  to  an  experimentally 
produced  fracture. 

The  utilization  of  this  data  and  its  correlation  to 
in  vivo  fracture  situations  has  been  discussed. 


•M 


VI . 


RECOMMENDATIONS 


This  section  briefly  describes  projected  efforts  that 
represent  logical  extensions  of  the  work  performed  in  this  con¬ 
tract.  The  curves  presented  in  Section  II  should  be  extended 
using  a  materials  testing  system  capable  of  loading  rates 
higher  than  those  obtainable  with  the  Instron.  It  would  also 
be  interesting  to  apply  a  greater  compressive  load  to  exagger¬ 
ate  the  effect  of  a  combined  axial  and  torsional  state  on  frac¬ 
ture  behavior.  Using  stress  analysis  techniques  a  model  could 
be  developed  to  predict  the  angle  of  fracture  and  to  explain 
why  a  more  complex  loading  state  produces  increased  magnitudes 
for  torque  and  energy.  Additional  specimens  should  be  tested 
to  better  define  the  role  of  cancellous  bone  in  energy 
absorption . 

The  efforts  described  in  Section  III  and  IV  should  be 
continued  and  expanded.  Specifically,  the  cross  sectional  and 
physical  properties  of  an  increased  number  of  human  and  animal 
bone  specimens  should  be  determined  to  examine  variations 
within  species.  Also,  a  pair  of  canine  radii  should  be  scanned 
to  determine  the  reason  for  consistant  fracture  at  the  proximal 
end  of  the  bone  as  described  in  Section  II.  The  relationship 
between  CT  number  and  density  must  be  derived  explicitly  and 
verified.  Correlation  of  in  vivo  and  in  vitro  fracture  sites 
with  cross  sectional  and  mechanical  properties  should  be 
con  tinned . 
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